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1.  Introduction. 

An  implicit  differential  algebraic  equation  (DAE)  has  the  form  i)  =  0  of  a  general 
ordinary  differential  equation  (ODE)  but  involves  a  (sufficiently  smooth)  mapping  F  : 
R”  X  R"  R”,  F  =  F{x,p).  for  which  rank  DpF{x,p)  =  r  <  n  is  constant  but  not  full 
in  F~^(0).  The  terminology  derives  from  the  fact  that,  in  practice,  the  simplest  DAE 
problems  are  of  the  form 

g(xi,X2)  =  0, 

i2  =  fixi,X2),  (a:i,X2)  €  R"“'' X  R’’, 

where  g  and  /  map  into  R"”"'  and  R*",  respectively,  and  thus  exhibit  clearly  a  differential 
and  an  algebraic  part. 

Because  of  the  rank  condition  a  DAE  does  not  reduce  locally  to  an  explicit  ODE  via  the 
implicit  function  theorem.  Nevertheless,  local  reduction  of  a  DAE  to  an  explicit  ODE  on  a 
sitbmanifold  of  R"  was  recently  shown  to  be  feasible  under  general  assumptions  ruling  out 
certain  geometric  singularities  ([RR91b]).  The  combination  of  this  reduction  procedure 
with  standard  ODE  results  yields  an  existence  and  uniqueness  theory  for  “nonsingular’’ 
DAE  s  ([RR91a],[RR91b])  which  is  applicable  to  many  problems. 

While  the  above  stresses  strong  analogies  between  DAE’s  and  classical  (that  is,  exphcit) 
ode’s,  important  differences  exist.  For  instance,  a  DAE  F{x,x)  =  C  need  not  have  a 
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solution  satisfying  the  initial  condition  x(0)  =  a;o,x(0)  =  po  vmder  the  sole  hypothesis  that 
F(xo,po)  =  0;  in  fact,  further  compatibility  conditions  are  usually  required  of  (xo,po)  for 
such  a  solution  to  exist.  Next,  solutions  of  DAE’s  may  exhibit  features  that  solutions  of 
explicit  ode’s  cannot  possess.  As  eui  example  with  n  =  2,  consider  the  DAE 

(  +  X2  =  0, 


(1.1) 


< 


with  the  initial  condition  x(0)  =  (1,-1)  (and  hence  x(0)  =  (—1/2,1)).  This  problem  has 
the  unique  solution  x(t)  =  ((1  —  —  1),  which  cannot  be  continued  beyond  t  =  1 

despite  the  fact  that  x(l)  =  (0,0)  exists  and  lim  x(t)  =  a'(l).  K,  instead,  x(t)  were 
characterized  as  a  solution  of  an  explicit  ODE  x  =  h(x)  with  continuous  h,  it  would  follow 
from  the  initial  condition  x(l)  =  (0,0)  that  x(^)  can  be  continued  beyond  t  =  1. 

The  point  (0,0)  where  the  solution  x(t)  =  ((1  —  —  1)  of  (1.1)  terminates  in  finite 

time  t  =  1  is  an  impassse  point,  as  they  are  called  in  the  engineering  literature,  especially 
in  connection  with  nonlinear  LRC  circuits  (see  e.g.  [C69],  [CD89]). 

We  concentrate  here  on  quasilinear  equations 


(1.2)  A(x)x  =  G(x), 

where  A  =  R”  £(R”)  satisfies  rank  A{x)  ~  r  <  n  whenever  G(x)  €  rge  -4(x).  This 
corresponds  to  the  choice  F(x,p)  =  A(x)p  —  G{x)  when  F  is  linear  in  p.  The  purpose 
of  this  paper  is  to  provide  a  mathematical  characterization  of  the  impasse  points  of  such 
quasilinear  DAE’s  which,  to  date,  is  lacking  in  general.  However,  the  recent  work  by 
Chua  and  Deng  [CD89]  contains  partial  results  for  the  special  case  when  A  is  a  constant 
projection. 

Despite  their  rather  special  structure,  quasilinear  DAE’s  are  relevant  in  a  large  number 
of  physical  problems.  Obviously,  t.he  simple  2x2  example  (1.1)  fits  into  this  framework. 
Moreover,  nonautonomous  problems  A{t,x)x  =  G{t,x)  can  also  be  brought  into  the  form 
(1.2)  by  adding  the  equation  t  =  1. 

Although  impasse  points  have  no  analogue  in  the  solutions  of  explicit  ODE's.  they  are 
strongly  reminiscent  of  “standard  singular  points”  of  singular  ODE’s.  Singular  ODE'.s 
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differ  from  DAE’s  in  that,  when  considered  in  the  form  (1.2),  the  matrix  A{x)  fails  to  have 
full  rank  n  only  at  exceptional  points.  As  shown  in  [RS9]  the  most  frequently  encountered 
type  of  singularity  for  a  singular  ODE  is  a  standard  singular  point  a:*  where  exactly  two 
solutions(^)  either  emanate  or  terminate.  This  is  exactly  the  case  in  the  simple  example 
(1.1)  where  t  i— »  (— (1  —  #—  1)  is  a  solution,  different  from  x{t)  =  ((1  — f  —  1),  that 

terminates  at  (0,0)  for  t  =  1.  This  strongly  suggests  that,  in  the  vicinity  of  an  impasse 
point,  a  DAE  (1.2)  can  be  reduced  to  a  singular  ODE  with  a  standard  singular  point. 
In  our  presentation  here,  the  possibility  of  performing  such  a  reduction  will  serve  as  an 
equivalent  mathematical  definition  for  impasse  points. 

For  convenience,  the  main  results  about  singular  points  are  summarized  in  Section  2. 
A  reduction  procedure  for  DAE’s,  which  is  essentially,  but  not  completely,  a  specialization 
of  the  method  in  [RR91b]  to  the  quasilinear  case,  is  described  in  Section  3.  The  value  of 
this  procedure  is  to  reduce  the  index  of  a  DAE,  and  in  particular  to  make  an  index  1  DAE 
into  an  implicit  ODE  that  may  or  may  not  be  singular.  The  short  Section  4  summarizes 
the  existence  theory  of  [RR91b]  for  the  index  1  quasilinear  DAE’s. 

Impasse  points  are  defined  in  Section  5.  We  have  chosen  to  give  a  coordinate-free  but 
rather  abstract  definition  based  upon  the  concept  of  intrinsic  derivatives  of  a  vector  bundle 
morphism.  This  definition  is  equivalent  to  a  simpler,  coordinate  dependent  one,  which  is 
useful  in  most  practical  applications  (Lemma  5.1). 

In  some  simple  but  frequent  cases,  impasse  points  turn  out  to  coincide  with  foldpoints 
of  a  specific  manifold  relative  to  some  splitting  of  the  ambient  space.  These  matters  are 
considered  in  Section  6  and  relate  to  the  results  of  [CD89]. 

Section  7  is  devoted  to  two  examples.  The  first  example  is  a  concrete  one  taken  from 
electrical  network  theory.  The  second  example  shows  that  one-parameter  stationary  prob¬ 
lems  can  be  reformulated  as  DAE’s  and  is  meant  to  emphasize  the  fact  that  the  study 
of  DAE’s,  including  impasse  points  and  singularities  (an  aspect  not  touched  upon  here) 
encompasses  all  the  issues  usually  investigated  in  the  framework  of  bifurcation  theory. 

The  theory  can  easily  be  extended  to  higher  index  problems  (provided  that  the  index 
is  well  defined).  Here,  the  validity  of  our  results  for  nonconstant  A(.t)  in  (1.2)  is  essential, 

f').\IoiJulo  tran.slations  in  tiinc  sinc^  (11)  is  autonomous. 
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even  if  -4(a:)  is  constant  in  the  given  higher  index  problem.  Some  related  comments  can 
be  found  in  Section  8. 

2.  Standard  Singular  Points  of  Singular  ODE’s:  A  Brief  Review. 

The  material  presented  here  is  taken  from  [R89]  but  our  notation  corresponds  to  that 
used  below  in  this  paper.  We  begin  by  introducing  the  following  terminology: 

Definition  2.1.  Consider  the  implicit  ODE 

(2.1)  = 

where  Aj  :  — >  £(]R'")  and  Gi  :  U'^  R''  are  mappings  defined  on  some  open  set 

C  A  point  ^  E  is  a  regular  point  of  (2.1 )  if  rank  Aj  (^)  =  r  and  a  singular  point 

if  rank  Ai((f )  <  r  hut  is  a  limit  point  of  regular  points  of  (2.1 ). 

Clearly,  for  any  regular  point  and  given  t*  the  initial  value  problem 

(2.2)  =  au)  =  u 

has  a  unique  solution  in  a  neighborhood  of  But  already  simple  examples  show  that  for 
a  singular  point  the  existence  and  number  of  soh'.tions  of  (2.2)  in  the  neighborhood  of 
if,  depends  strongly  on  the  character  of  the  singularity  at  that  point.  The  most  frequent 
and  also  simplest  case  is  the  following; 

Definition  2.1.  A  singular  point  f  G  U'"  of  (2.1)  is  a  standard  singular  point  if 

(2.3)  dimker  Ai(f )  =  1  (  i.e.  rank  Ai(f )  =  r  —  1), 

(2.4)  G,(f)^  rgeAAf). 

(2.5)  (DAi(f)7/)7/ ^  rgeAi(f),  Vr/ G  ker  Ai (f )  x  {0}. 

In  particular,  if  (2.1)  has  a  standard  singular  point  then  such  points  will  exist  for 
every  problem  obtained  by  a  sufficiently  small  perturbation  of  (2.1).  Moreover,  the  set  of 
standard  singular  points  of  (2.1)  is  either  empty  or  forms  a  hypersurface  of  and.  at 


least  in  the  case,  the  standard  singularities  are  the  only  singularities  that  satisfy  the 
latter  condition  for  generic  choices  of  and  Gi  (see  [R89]). 

Remark  2.1:  The  fact  that  standard  singular  points  of  (2.1)  lie  on  a  hypersurface, 
if  they  exist  at  all,  has  an  important  consequence  regarding  the  dynEimics.  Indeed,  if  a 
solution  of  the  initial  value  problem  (2.2)  for  some  regular  point  tends,  at  some  later 
time,  to  a  standard  singular  point,  then  the  same  will  happen  also  with  any  trajectory 
passing  through  a  point  sufficiently  close  to  ^o-  In  other  words,  standard  singular  points 
cannot  be  by-passed  by  simply  perturbing  the  initial  condition.  □ 

If  is  a  standard  singular  point  of  (2.1),  the  condition  (2.4)  makes  it  impossible  for 
any  function  ^  -|-  T)  — >  R’'  with  ^(f*)  =  to  be  a  solution  of  (2.1)  since  the 

relation  )^(t,)  =  )  cannot  hold.  This  remark  suggests  the  following  definition. 

Definition  2.3.  Let  be  a  standard  singular  point  of  (2.1).  A  solution  of  the  initial 
value  problem  (2.2)  for  given  t,  €  R  is  a  continuous  function  ^  :  J  C  R^  — >  b’’'’  on  either 
J  =  [t*,  t*  +  T)  or  J  =  (t,  —  T,  f*]  for  some  T  >  0  which  is  of  class  on  =  J  } 
and  satisfies  ^(t.)  =  and  Ai{^{t))^(t)  =  Gi(^(t))  fort  6 

Evidently  the  two  conditions  (2.4)  and  (2.5)  are  equivalent  with 

{(DAi{OT))T},fj)  ^  0, 

(2.6)  V77  €  kerAif^)  \  {0},  rj  e  ker  Ai(0^  \  {0} 

where  (•,  •)  denotes  the  natural  inner  product  of  R'^.  Since  the  form  fj)  is  continuous 

and  quadratic  in  rj  as  well  as  fj  it  must  have  either  a  positive  or  negative  value  for  all  pairs 
of  nonzero  vectors  (7,77)  in  ker^i(^)  x  ker.4i(^)^  if  only  this  holds  for  one  such  pair. 

The  main  result  ([R89,  Theorem  5.1])  on  the  existence  of  solutions  near  standard  sin¬ 
gular  points  can  then  be  phrased  as  follows; 

Theorem  2.1.  Let  be  a  standard  singular  point  of  (2.1 ).  Then,  given  G  R.  the  initial 
value  problem  (2.2)  has  exactly  two  solutions  which  are  both  defined  on  J  =  -|-  T) 

or  on  .7  =  (f.  —  r,  f*|  for  some  T  >  0  if  a(^,)(r/,  ^)  >  0  or  a(i^,)(77, 77)  <  0  for  some  pairs 
of  nonzero  vectors  (r/,ij)  G  ker.4i(if*)  x  ker>li(^»)^,  respectively.  Furthermore.  ||(f(f)|| 
tends  to  oc,  as  f  €  -7  \  {<*}  tends  to  f,. 
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It  follows  from  Theorem  2.1  that  a  solution  of  (2.1)  emanating  from  some  regular  point 
can  reach  a  standard  singular  point  at  some  later  time  only  if  the  quantity  Q;(^*)  of  (2.6) 
is  negative.  Standard  singular  points  (f*  with  positive  Q;(^»)  obviously  cannot  be  reached 
in  increasing  time.  Thus,  it  is  appropriate  to  introduce  the  following  terminology: 

Definition  2.4.  The  standard  singular  point  of  (2.1)  is  accessible  or  inaccessible  if 
if)  >  0  or  a(^t)(rj,  if)  <  0,  respectively,  for  some  pair  of  nonzero  vectors  tj,  fj  in 
ker  Ai(^*)  x  ker  Ai(^,)^. 

In  [RS9],  accessible  and  inaccessible  points  are  called  attracting  and  repelling,  respec¬ 
tively,  but  the  terminology  used  here  appears  to  be  more  appropriate.  Of  course,  by 
reversing  the  evolution  in  time,  inaccessible  points  become  accessible  and  vice-versa. 

It  follows  from  Theorem  2.1  that  accessible  standard  singular  points  are  reached  in 
finite  time  by  trajectories  emanating  elsewhere  in  U^.  Since  these  trajectories  cannot  be 
continuously  extended  beyond  these  points,  they  represent  “catastrophes”  for  the  solutions 
of  (2.1).  This  is  further  emphasized  by  the  fact  that,  as  mentioned  in  Remark  2.1,  no  small 
perturbation  of  the  initial  condition  (and/or  of  Ai  or  Gi)  will  affect  the  eventual  encounter 
of  such  points.  In  this  respect,  note  that  the  sign  condition  in  (2.6)  is  unchanged  if  Ai, 
G\  and  are  replaced  by  sufficiently  small  approximations. 

Although  they  cannot  be  reached,  inaccessible  standard  singular  points  may  also  have 
drastic  effects  on  the  dynamics  but  not  in  the  form  of  catastrophes  (see  again  [R89]). 

3.  A  Reduction  Procedure  for  Differential-Algebraic  Equations. 

In  this  section,  we  consider  DAE’s  of  the  form 

(3.1)  A{x)i  =  G{x),  r  €  K" 
with  the  following  properties: 

Assumption  3.1.  For  some  open  subset  U”  C  R"  the  mappings  A  :  — >  £(R")  and 

G  :  U"  — >  R"  are  of  class  on  and  for  some  fixed  integer  0  <  r  <  n,  \vc  have 

(3.2)  {.r  G  t"".  G{x)  G  rgc  A(.r)}  =>  rank  A(x)  =  r 


tt 


and  the  mapping 


(3.3)  (x,p)  6  LI"  X  R"  ^  A(a:)p  -  G(x)  €  R”, 
is  a  submersion. 

The  submersion  assumption  for  (3.3)  implies  that  the  mapping 

(3.4)  {h,q)  €  R"  X  R”  I — >  {DAix)h)p  +  A{x)q  -  DG(x)h  e  R" 
is  surjective  for  every  (x,p)  €  x  R".  Hence  the  set 

(3.5)  M  =  {(x,p)  e  U’^  X  R"  :  A{x)p  —  G(x)  =  0} 

is  a  closed  n-dimensional  submanifold  of  x  R". 

Our  aim  is  to  show  that  (3.1)  can  be  reduced  to  the  similar  form 

(3.6)  Adi)i  =  Gi{0. 

where  now  Ai  and  G\  take  values  in  £(R'')  and  R'",  respectively.  In  essence,  this  reduction 
represents  a  simplification  in  the  quasilinear  case  (3.1)  of  the  theory  developed  in  [RRQlb] 
for  general  implicit  DAE’s  F(x,x)  =  0.  However,  there  is  a  special  feature  which  will  be 
of  crucial  importance  for  the  later  discussion.  While  the  reduction  of  [RRQlb]  is  valid  only 
locally  in  the  vicinity  of  a  point  (x*,p, )  €  i^~V0),  the  reduction  here  is  local  only  in  the 
first  variable.  This  is  due  to  the  linearity  of  (3.1)  with  respect  to  the  derivative  x  and, 
roughly  speaking,  will  allow  us  to  analyze  phenomena  involving  “infinite”  p*  as  must  be 
done  to  discuss  impasse  points. 

We  consider  the  set 

(3.7)  ir  =  {j  e  b'”  :  G{x)  €  rge  A(x)), 

and  note  that  (x,p)  G  M  for  some  p  G  R"  if  and  only  if  x  G  IT  and  hence 

(3.8)  IT  =  7r{M), 

where  "  :  f’"  x  R"  — ♦  is  the  projection  onto  the  first  factor.  Then  we  have; 


Proposition  3.1.  The  set  W  is  an  r- dimensional  submanifold  ofU'^.  In  addition,  if 

the  set  {x  €  '■  rankA(x)  =  r}  is  closed  in  U"  then  W  is  closed  in  Z7”. 

Proof:  By  Assumption  3.1,  rank  A(x)  =  r  is  independent  of  x  €  W  which  implies  that 
rank  is  constant  and  equed  to  r.  Indeed,  for  (x,p)  €  M,  the  mapping  T(j.  p)(7r|j^^ )  is 
the  restriction  of  tt  to  T(j.^p)M  and,  since  T(j.p)M  is  the  null-space  of  the  mapping  (3.4), 
we  have 

kerr(,,p)(7r|^)  =  {(0,g)  €  R"  :  A{x)q  =  0}  =  {0}  x  ker  A(x). 

Thus,  for  (x,p)  €  M,  and  hence  x  €  W,  we  obtain  dimker  T(j.  p)(7r|j^^ )  =  n  —  r  and  therefore 
rank  =  r  as  desired. 

In  view  of  this  result  and  the  subimmersion  theorem  (see  e.g.  [D70],  [AMRSS])  it  suffices 
to  show'  that  :  M  ^  W  is  open  to  prove  that  W  is  an  r  -  dimensional  submanifold 
of  U^.  Let  Mo  be  any  open  subset  of  M.  In  order  to  prove  that  7r(Mo)  is  open  in  W 
let  xo  €  7r(Mo),  so  that  there  is  some  po  €  R"  with  (xo,po)  G  Mo-  We  construct  a 
continuous  function  /  :  Wo  —*  M  on  some  open  neighborhood  Wo  C  W  of  xo  such  that 
/(xo)  =  (io,po)  and  tt  o  /  =  idwo-  The  continuity  of  /  and  (xo,po)  G  Mo  then  imply 
that  f~^(Mo)  is  an  open  neighborhood  of  xo  in  Wo  and  hence  in  W.  On  the  other  hand, 
TT  0  /  =  idw'o  implies  C  7r(Mo).  Thus  it  follows  that  7r(Mo)  is  a  neighborhood  of 

Xo  and  therefore  that  tt{Mo)  is  open  since  xo  wa^  arbitrary. 

Before  constructing  /  let  us  observe  that,  in  general,  if 

(3.9)  F  =  {x  €  :  rank  A{x)  =  r) 

and  Xo  €  V  is  a  given  point,  then  there  exist  an  open  neighborhood  Vq  of  in  V  and 
continuous  functions  e,  :  Vo  —*  R"  such  that  {ei(x),---  ,er(x)}  is  an  orthonormal  basis 
of  rge  A{x),  Vx  e  Fq.  Indeed,  let  uoi,---  -Uor  be  any  basis  of  rge  .4(xo).  Then  there 
are  linearly  independent  vectors  Wi,---  ,Wr  €  R'’  for  which  uoi  =  -4(xo)te,,l  <  i  <  r. 
A  straightforward  contradiction  argument  now  shows  that  the  vectors  t/,(x)  =  A{x)u',, 
i  =  1.  •  •  •  .r.  remain  linearh'  independent  for  all  x  sufficiently  close  to  .To  in  Thus,  for 
all  x  in  some  small  open  neighborhood  I'o  of  xo  in  F  the  vectors  tii(x),---  ,»,(x)  form 
a  basis  of  rge  A{x}  .since  rank  A(x)  =  r.  Now',  by  applying  the  Gram-Schmidt  process 
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to  Ui{x),---  ,iir(x),  we  obtain  an  orthonormal  basis  ,er(x)}.  The  functions 

Ci  :  Vo  R"  are  continuous  since  the  functions  x  it,(ar)  are  continuous  and  the  Gram- 
Schmidt  process  involves  only  continuous  operations.  For  future  use  we  note  that  the 
orthogonal  projections  P(x)  €  £(R")  onto  rge  .A(x)  depend  continuously  upon  x  €  TFo 
since  P{x)  can  be  expressed  as  a  sum  of  dyadic  products  of  the  vectors  ei(x),  •  •  •  ,  er(x). 

Clearly,  the  above  arguments  can  be  repeated  verbatim  to  obtain  continuous  functions 
e;  :  ^  z  =  1,  •  •  ■  ,  r,  such  that  {ei(x),  •  •  •  er(x)}  is  an  orthonormal  basis  of  rge  A{x)'^ 

(after  shrinking  Vo  if  necessciry).  In  particular,  the  orthogonal  projection  P{x)  G  £(R") 
onto  rge  -4(x)^  depends  continuously  upon  x  G  Vq- 

We  return  to  the  specific  case  when  xo  €  W.  With  the  constructed  neighborhood  To 
we  obtain  the  open  neighborhood  TT^o  =  W  H  T-'o  of  xo  in  TV.  For  x  G  TVo,  the  equation 
.4(x)p  =  G{x),  p  G  {ker^(x)}-‘-  =  rge  -4(x)^,  has  a  unique  solution  p  =  p{x).  The 
components  p,(x)  of  p(x)  in  the  basis  ei(x),  1  <  z  <  r,  satisfy  the  (invertible)  linear 

r 

system  0'ij(x)pj(x)  =  7j(x),  1  <  z  <  r,  where  aij(x)  =  (e,(x), >l(x)e_,(x))  and 

j=\ 

lx{x)  =  {G(x),ei(x))  are  continuous  on  TVo-  Therefore  the  components  pi,  and  hence  p 
itself,  are  continuous  functions  on  TT4). 

We  construct  the  desired  function  /  in  the  form  /(x)  =  (x,p(x)  +  g(x))  where  q  :  TT'o  — > 
R"  is  continuous  and  satisfies  ^(x)  G  kerA(x)  for  x  G  TVq,  and  q{xo)  =  po  -p(xo)  € 
ker^(xo).  Clearly  then  /(TVq)  C  M,  /(xq)  =  (xo,po)  and  n  o  f  —  idwo-  In  order  to  find 
q  it  sufiices  to  note  that,  in  view  of  the  previous  remarks  Q  —  I  —  P  :  Wq  — >  £(R")  is 
continuous  and  Qix)  is  the  orthogonal  projection  onto  ker.4(x)  =  {  rge  /l(x)^}-‘-.  Thus, 
q{x)  —  Q(x)(po  —  p(xo))  satisfies  the  desired  conditions. 

Fo  the  completion  of  the  proof,  we  must  show  that  TV  is  closed  in  P"  whenever  the 
set  V  of  (3.9)  is  closed  in  C”.  Let  {x^}it>i  be  a  sequence  in  TV  with  lim  Xk  =  xq  G  C". 

fc— >oo 

Then  we  have  xo  G  V  since  TT^  C  V  and  V  is  closed  in  C/".  Let  To  C  T'  be  an  open 
neighborhood  of  xq  in  V  such  that  P  :  Tq  — >  £(M")  is  continuous  where  P(x)  denotes  the 
orthogonal  projection  onto  rge  A(xy,  the  existence  of  such  a  neighborhood  Vo  was  shown 
above.  Because  of  Xk  6  TV  we  have  P{xk)G{xk)  =  G(xk)  for  k  >  1.  In  the  limit  as 
V  — >  oc  it  then  follows  that  P(xo)G(xo)  =  G(xo)  which  means  that  G(xo)  G  rge  .4(xo) 
and  therefore  xo  G  IT'.  □ 
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Remark  3.1:  Proposition  3.1  has  no  global  analog  for  general  implicit  DAE’s  F(x,  x)  = 
0,  although  a  local  version  exists  (see  [RR91b]).  □ 

We  are  now  in  a  position  to  describe  our  reduction  procedure.  As  a  motivation,  note 
that  if  J  C  K  is  an  open  interval  and  x  :  J  — »  [7"  a  solution  of  (3.1)  then  necessarily 
x{t)  G  W,  yt  G  J.  Since  W  is  a  manifold  this  implies  that  {x{t),x{t))  G  TW,  Vt  ^  J  where 
here  and  subsequently  TW  is  viewed  as  a  subset  of  TR”  =  R”  x  R".  Since  (x(t),  x{t))  G  M 
it  follows  that  (x(i),x(t))  G  TW  D  M,  Vt  G  J-  The  reduction  of  (3.1)  to  the  form  (3.2) 
now  is  an  easy  consequence  of  the  characterization  of  TW  Pi  M  given  below. 

Let  X,  G  W.  By  Proposition  3.1  there  exist  open  subsets  O  C  W  and  U''  C  R'^  and  a 
mapping  ip  :  R'"  which  is  a  diffeomorphism  of  onto  O.  The  tangent  bundle 

TO  is  the  image  of  the  mapping 

(3.10)  (VP(0,  DifiOr})- 
From  TO  =  {(x.p)  G  TW  :  x  G  O}  we  obtain 

(3.11)  {(x,p)GrwnM:  xGO}  =  rc»nM, 

and,  using  (3.11)  and  the  fact  that  TO  is  the  image  of  the  mapping  (3.10),  that 

(  x  =  v?(0,  P  =  DifiOv, 

(3.12)  {(x,p)  G  TW  n  M,  xeO}^  < 

[  A{ip{0)Dip{0v  -  (?(<^(e))  =  0, 

for  some  pair  (^,77)  G  x  R''. 

By  definition  of  ip  and  W ,  G(i,5(^))  G  rge  A((7?(^)),  G  which  implies  that 

(3.13)  -  G{ip{0)  €  rge  A(vp(0),  G  GF 

As  we  shall  see  now.  (3.13)  allows  us  to  replace  the  equation  in  R"  appearing  on  the 
right-hand  side  of  (3.12)  by  a  similar  equation  in  R'". 

Let  Z,  be  a  given  complement  of  rge  A(x,).  Then  Z,  is  also  a  complement  of  rge  A(x) 
for  all  X  G  W  sufficiently  close  to  x,.  This  can  be  shown,  for  instance,  by  using  the  fact 
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that  the  orthogonal  projection  P{x)  onto  rge  A[x)  depends  continuously  on  x  €  W ,  as 
we  saw  in  the  proof  of  Proposition  3.1.  Now  suppose  that  there  is  a  sequence  {xjt}ii.>i  in 
TP  tending  to  x,  such  that  Z,  is  not  a  complement  of  rge  yl(xt).  Since  rank  A{xk)  —  r 
and  dimZ,  =  n  —  r,  there  exists  a  unit  vector  uk  €  rge  A{xk)  D  Z,  and,  by  extracting  a 
subsequence,  we  may  assume  that  lim  Uk  =  u.  Clearly,  u  is  a  unit  vector  on  Z,  and,  since 

k—^oo 

P{xk)uk  =  Uk  it  follows  that  P(x,,)u  =  u;  that  is,  u  6  rge  A(x,)  which  is  absurd. 

Thus,  by  shrinking  O  and  with  it  f/'’,  if  necessary,  we  can  ensure  that  Z,  is  a  complement 
of  rge  .4((/?(if)),  for  ^  €  17’".  Equivalently,  the  projection  operator  P»  onto  rge  -4(x»)  relative 
to  the  decomposition 

(3.14)  R"  =  rge  A(x.)0Z. 

is  a  linear  isomorphism  of  rge  .4(<,5(^))  onto  rge  A{x^)  for  all  ^  e  17’'.  By  (3.13)  this  implies 
that  for  any  pair  (^,7?)  €  IP'  '<  R’’  we  have  A{(p{^))D(i:>{Ov  -  G{(p{^))  =  0  if  and  only  if 
P»A{!Tp{())Dip{^)r)  —  P,G((^(0)  =  0.  Using  this  equivalence  in  (3.12),  we  infer  that 

(  X  =  P  = 

(3.15)  {{x,p)  eTWnM.  X  eO}  ^  < 

U4i(Ori-Gi(O  =  0, 

for  some  pair  (^,r/)  €  17'^  x  R’’,  where  we  have  set 

(3.16)  A,{0  =  P,A{^{0)D^{0 

(3.17)  =  PM^'iO)- 

Note  that  by  identifying  rge  i4(x.)  with  R’’  we  see  that  the  mappings  Ai  and  Gi  map  into 
£(R’')  and  R’",  respectively,  and  are  of  class  CU 

Suppose  now  that  .7  C  R  is  an  open  interval  and  x  :  J  — >  7’"  a  C*  solution  of  (3.1).  As 
noted  earlier,  we  have  x(J)  C  TP.  If  it  so  happens  that  x(.7)  C  G  so  that  (x(f).x(/))  6 
TIP  n  M  and  x{t)  €  O  for  all  f  £  J.  then  it  follows  from  (3.15)  that  the  C*  function 
f  :  .7  — >  7’’'  defined  by  if(f)  =  o  x(t)  satisfies 


(3.1S) 


-4,(^(i))ao  =  c;,(e(/)).  vf  €  J. 


with  -4i,Gi  given  by  (3.16)  and  (3.17),  respectively.  Conversely,  for  everj^  solution 
(^  :  J  — »  L'”'’  of  (3.10),  x{t)  =  is  a  solution  of  (3.1).  Henceforth,  Equation  (3.18) 

will  be  referred  to  as  the  reduction  of  (3.1)  near  x*. 

4.  Existence  Theory  for  Nonsingular  DAE’s  of  Index  1. 

The  results  presented  in  this  section  are  particular  cases  of  the  general  theory  in  [RROlb] 
but  the  main  existence  and  uniqueness  Theorem  4.1  can  here  be  given  an  independent  and 
much  simplified  treatment  by  talcing  advantage  of  the  special  quasilinear  form  (2.1)  of  the 
DAE  and  by  confining  attention  to  index  1  problems  as  defined  below.  We  retain  the 
hypotheses  and  notation  of  Section  3. 

Definition  4.1.  The  quasilinear  system  (3.1)  (satisfying  Assumption  3.1)  is  a  nonsingular 
DAE  of  index  1  if 

(4.1)  (x  €  C",  G(x)  e  rge  A(x)|j.^,y  }(^)  ranic  i4(x)|j,^^  =  rank  A(x)(=  r). 

For  nonsingular  DAE’s  of  index  1  the  local  existence  and  uniqueness  theory  is  very 
similar  to  that  for  ODE’s: 

Theorem  4.1.  Let  (3.1)  be  a  nonsingular  DAE  of  index  1.  Then,  for  any  given  x,  6 
IF]  =  7r(TH'  n  A/)  C  W  and  t*  G  K,  there  exists  an  open  interval  J  containing  t.  and  a 
unique  solution  x  :  J  —y  U'^  of 

(4.2)  A{x)x  =  G{x),  x(t.)  =  x.. 

Moreover,  no  solution  of  (4.2)  exists  for  x,  ^  TFi . 

Proof:  In  Section  3  it  was  already  noted  that  when  x  ;  J  — >  U"  is  a  solution  of  (3.1) 
tlu'ii  x(t))  €  TW  n  M.  Vt  £  J.  and  hence  x(.7)  C  W\  —  7r(ril’  fl  A/).  In  particular. 

(4.2)  cannot  have  a  C*  solution  if  .r,  ^  Hq . 

Suppose  now  that  .r*  £  IF]  C  IF  so  that  (x*,p*)  £  TW  fl  A/  for  some  p*  £  M".  Then 
p,  £  Tj.,  IF  and  .4(.r*)p*  =  G'(.r,)  imidies  that  G'(,r,)  £  rge  .4(.r»)|j,^  and  therefore,  by 

(-'jlloiKo  X  G  M';  SCO  (3.7). 
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(4.1),  that  rank  w  =  Equivalently,  this  shows  that 

(4.3)  A{x^)  e  GL{  rge  A(a:*)  ). 

The  reduction  procediure  of  the  previous  section  ensured  the  existence  of  a  diffeo- 
morphism  9  of  an  open  subset  of  R'’  onto  an  open  neighborhood  O  oi  in  W  such 
that  on  O  (4.2)  is  equivalent  to 

(4.4)  = 

where  1^,  =  <^“^(ar*),  Ai{()  =  PtA(<f{^))D<f(^),  Gi{^)  =  P.G((,?(^))  and  P*  denotes  a 
projection  operator  onto  rge  ^(a;,). 

For  ^  =  1^.  we  have  ^i(^,)  =  A(x^)D(p{^^)  6  GZi(R'",  rge  A{x^))  by  (4.3)  and 

Pi,^(C*)  €  GT(R'', Tj;,  W).  Thus,  in  the  vicinity  of  (4.4)  is  equivalent  to  the  explicit 
ODE  ^  =  ^i(^)“^Gi(^),  ^{t»)  =  ^*.  Hence  the  classical  ODE  theory  guarantees  that  on 

some  open  interval  J  containing  t,  there  exists  a  imique  solution  (^  :  J  — >  P’’  of  (4.4). 
Clearly  then  a:  =  (^  o  is  the  unique  solution  of  (4.2).  □ 

Remark  4.1:  It  is  easily  checked  from  the  above  proof  that  the  solution  x  in  Theorem 
4.1  is  actually  of  class  C^.  □ 

The  Definition  4.1  of  a  nonsingular  DAE  of  index  1  does  not  rule  out  the  existence 
of  points  ar*  €  W  where  rank  A(x*)|j.^  ^  <  r.  Such  points  do  not  belong  to  the  set 
ITj  =  Tr(TW  n  M)  and  hence,  by  Theorem  4.1,  no  solution  to  the  corresponding  initial 
value  problem  (4.2)  exists.  Nevertheless,  on  some  interval  J  =  —  T.  t,]  or  J  =  [f,,<*  +  T] 
with  T  >  0,  continuous  functions  may  well  exist  which  are  differentiable  on  J  \  {t,}  and 
satisfy  there  A{x{t))x{t)  =  G(x(t)).  The  impasse  points  discussed  in  the  next  section  allow 
for  solutions  of  this  type. 

5.  Impasse  Points  of  Nonsingular  DAE’s  with  Index  1. 

First,  we  recall  the  definition  and  basic  properties  of  the  intrinsic  derivative  of  a  vector 
bundle  morphism  (.see  e.g.  [AGVS5],  [GG73]).  Let  A'’  be  a  manifold  and  E.F  vector 
bundles  with  base  A'’.  A  vector  bundle  morphism  p  :  E  F  can  always  be  viewed  as  a 
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section  p  :  X  — >  Horn  {E,F)  of  the  vector  bundle  Horn  {E,F)  with  base  X.  Thus,  for 
X  E  X ,  we  have  p{x)  =  {x,  B(x))  where  B{x)  €  C(Ej.,  Fx)  and,  as  usual,  E^  and  F^  denote 
the  fibers  of  E  and  F,  respectively,  above  x  E  X.  If  p  is  differentiable,  T^p  is  a  linear 
mapping  TxP  :  T^X  — +  Tp(x)  Horn  (E,F).  As  Tp^^)  Horn  {E,F)  ~  T^X  x  C{Ei,Fx)  and 
p  is  a  section,  we  have  {Txp)h  =  (h,  B{x)h)  E  T^X  x  C{Ex,  V/i  E  T^X,  where  then 
B{x)ECiTxX,C{Ex,Fx)). 

The  intrinsic  derivative  of  p  at  a:  €  X  is  defined  as  the  mapping  (loc.  cit.) 

(5.1)  h  E  TxX  I — >  E  £(kerB(ar),  coker  B{x)), 

where  cr(a:)  :  coker  B(x)  —  F^/rge  B{x)  is  the  canonical  projection. 

When  X’  is  an  open  subset  of  R'’  and  F  =  X  x  R",  F  =  X  x  R*"  then  J3  is  a  mapping 
B  :  X  — +  £(R”,R'”)  and  B(x)  is  just  the  derivative  DB{x).  Hence  the  intrinsic  derivative 

(5.1)  becomes 

(5.2)  h  eW  I — >  <7(a')[FF(x)/i]|j^^^g^^^  €  £(ker  J5(x),  coker  B{x)), 

where  a{x)  :  R*"  — >  coker  B(x)  =  R'"/rge  B(x)  is  the  canonical  projection.  It  is  well 
known  that,  in  general,  the  calculation  of  the  intrinsic  derivative  (5.1)  reduces  to  the 
special  case  (5.2)  via  local  trivializations  of  the  vector  bundles  E  and  F  (loc.  cit.). 

In  our  setting,  a  useful  choice  is  given  by  X  =  W,  E  =  TW,  F  =  R\,v(A)  where 

(5.3)  Rw(A)  =  {(x,9)  eW  X  R";  q  E  rge  A(x)}. 

As  rank  A(x)  =  r  is  independent  of  x  e  W,  it  is  indeed  easily  seen  that  R\v{A)  is  a  vector 
bundle  with  base  W  and  r-dimensional  fiber  rge  A(x),  x  E  W.  We  set 

(5.4)  A\v  :  (x,p)  E  TW  * — >  A\y{x,p)  —  (x,  A(x)p)  E  Fu’(A), 
where,  alternatively.  Aw  may  be  viewed  as  the  section 

(5.4')  Aw  :  X  E  H'  i — >  (x.  A(.t)|^^ )  E  Horn  (TW,  R\v{A)). 
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Because  of  the  embedding  rge  A{x)  C  K",  coker  =  rge  A{x)l  rge 

may  be  identified  with 

(5.5)  coker  =  rge  ^(x)  D  (rge 

so  that  o{x)  :  rge  .4(x)  — >  coker  A{x)^J.^^^,  corresponds  to  the  restriction  to  rge  A{x)  of 
the  orthogonal  projection  onto  coker  A{x)\j.^^^ . 

Definition  5.1.  Let  (3.1)  be  a  nonsingular  DAE  with  index  1.  The  point  x,  G  is  an 
impasse  point  of  (3.1 )  if  the  two  conditions 

(5.6)  dimker  yl(x*)|^^  »v  =  1  (^■^-  rank  A(x,)ij.^^  =  r  —  1). 


(5.7) 


hold,  where  ' DA\v{x»)  denotes  the  intrinsic  derivative  of  the  bundle  morphism  A\v  at  x*. 
The  impasse  point  x*  is  said  to  be  accessible  or  inaccessible  if 


(5.S) 


a(^)(e,e)  =  (G(x.),  e)(('£>i4n/(x.)e)e,  e) 


is  <  0  or  >  0,  respectively,  for  some  pair  (e,e)  of  nonzero  vectors  in  ker.4(x,)|y^  x 
coker  j4(x»  )|^^  where  (-,■)  denotes  the  usual  inner  product  ofR". 

A  few  comments  are  in  order.  First,  as  in  the  case  of  (2.6)  it  follows  that  the  form 

(5.8)  must  have  either  a  positive  or  negative  value  for  all  pairs  of  nonzero  vectors  (e,e) 

in  ker  .4(.t,)|j-^  x  coker  .4(x,)|j.^  if  only  this  holds  for  one  such  pair  of  vectors.  Sec¬ 
ondly.  it  follows  from  (5.6)  that  'DA\v{x„)  is  a  linear  mapping  from  Ti^W  into  the  one- 
dinicnsional  space  C{ker  A{x„)\.j.^  ,  coker  .4(x,)|j.^  „. ).  Hence,  condition  (5.7)  is  equiva¬ 

lent  to  the  assumption  that 

(5.9)  ' DA\v(x»)e  ^  Q, 
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for  some,  or  equivalently  every",  nonzero  vector  e  €  ker  In  turn,  (5.9)  holds  if 

and  only  if  for  some,  or  equivalently  every,  nonzero  vector  e  G  coker  we  have 

(using  the  identification  in  (5.5)) 

(5.10)  (('DAu/(x*)e)e, e)  ^  0. 

On  the  other  hand,  because  of  x,  G  W  it  follows  that  G(x,)  G  rge  A(x.)  and  there¬ 
fore  by  (5.6)  and  the  hypothesis  that  the  DAE  (3.1)  is  nonsingular  (Definition  4.1)  that 
(G(x,),e)  7^  0  for  e  G  coker  >l(x*)|j.^  \  {0}.  Together  with  (5.10)  this  shows  that  the 

left-hand  side  of  (5.8)  is  nonzero.  Moreover,  it  is  obvious  that  its  sign  is  independent  of 
the  nonzero  vectors  e,  e. 

Remark  5.1:  It  should  be  emphasized  that  the  intrinsic  derivative  'DAw  cannot  be  re¬ 
placed  by"  the  much  simpler  DA  in  (5.7)  or  (5.8),  although  the  conditions  obtained  through 
the  substitution  formally  make  sense  but  are  nevertheless  irrelevant.  To  corroborate  the 
above  statement,  it  suffices  to  note  that  when  'DAw  is  replaced  by  DA  in  (5.7)  or  (5.8), 
these  modified  conditions  never  hold  when  A(x)  is  independent  of  x,  whereas  constancy"  of 
A{x)  is  not  an  obstacle  to  the  existence  of  impasse  points  as  we  will  see  in  the  next  section. 

The  following  lemma  relates  Definition  5.1  to  standard  singular  points  of  ODE’s  (see 
Section  2)  via  the  reduction  procedure  of  Section  3. 

Lemma  5.1.  Let  (3.1)  be  a  nonsingular  DAE  with  index  1.  The  point  x*  G  W  is  an 
accessible  or  inaccessible  impasse  point  of  (3.1)  if  and  only  if  =  ip~^(x»)  is  an  accessible 
or  inaccessible  standard  singular  point,  respectively,  of  the  reduction  (3.18)  of  (3.1)  locally 
near  x,. 

Note:  From  this  characterization,  it  follows  that  a  nonsingular  DAE  may  reduce  to  a 
singular  ODE.  Of  course,  this  is  not  due  to  any  ambiguity  in  the  choice  of  terminology  but, 
rather,  to  the  fact  that  the  reduction  procedure  may,  of  necessity",  give  rise  to  singularities. 

Proof:  The  proof  follows  from  the  calculation  of  the  intrinsic  derivative  'D.4u  (.t*) 
using  appropriate  trivializations  of  the  bundles  TW  and  /2iv(.4)  near  x,.  A  trivialization 
of  TW  near  x,  is  given  by  the  inverse  of  the  mapping  {f,T])  G  U'^  x  R’'  i — >  (9(0- 
where  V  is  an  open  subset  of  R*"  and  ^  a  diffeomorphism  of  onto  an  oi)en  neighborhood 
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C?  of  X,  in  W.  On  the  other  hand,  with  a  projection  operator  P*  onto  rge  A(i,),  the 
mapping  {x,q)  G  Rwi^.)  i — >■  P,?)  €:  U’"  x  rge  A(x,)  provides  a  trivialization 

of  Pn/(i4)  near  x*  after  shrinking  if  necessary.  Indeed,  as  noticed  in  the  proof  of 

Proposition  3.1,  we  have  P*  G  G£(rge  A(x),  rge  j4(x*))  for  x  G  W  sufficiently  close  to 

X.. 

The  local  expression  of  the  bundle  morphism  Aiy  in  (5.4)  becomes 

(5.11)  itri)  G  X  R’-  ^  {^,P,AMi))DipiOv)  e  x  rge  ^(x.)  ~  IP  x  R". 

As  mentioned  at  the  beginning  of  this  section,  the  corresponding  local  expression  of 
' DA\v{x^)  is  the  intrinsic  derivative  at  =  ip~^{x^)  of  the  morphism  (5.11)  of  trivial 
bundles  over  an  open  subset  of  R'".  Thus,  in  line  with  (5.2),  it  is  given  by  G  R'^  ' — > 

where 


(5.12)  A^U)  =  P,A{^{0)Dq:i0 

and  (t(^*)  ;  R'"  — »  R''/rge  Ai(^*)  is  the  canonical  projection. 

Conditions  (5.6)  and  (5.7)  amount  to 

(5.13)  dimker  Ai(^*)  =  1  (i.e.  rank  Ai(^,)  =  r  —  1) 

and  to  h  G  ker.4i(^, )  ' — >  /»,({.)  7^  respectively.  But.  in  view  of 

(5.13) .  it  is  clear  that  the  latter  relation  reads 

(5.14)  (DAi{^^)ri)T]  ^  rge  Ai(^,),  V77  G  ker  Ai(^,)  \  (0). 


Finally,  with 

(5.15)  6h(0  =  P.CTV(0).  VCGC^ 
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we  have  =  P*G(i*)  =  G(x,)  since  P*  =  /  on  rge  A(x,)  and  6  T-V’.  As 

rgeA(a:,)|j,^  =  rge  Ai(^,)  is  obvious  from  (5.12)  and  Dip{^^)  €  Tj.,  11^),  it 

follows  that 

(5.16)  G{xr,)  ^  rge  A(a;.)|j.^^j^,  ^  Gi(^.)  ^  rge  Ai(^,). 

Suppose  now  that  (3.1)  is  a  nonsingular  DAE  and  x,  6  is  an  impasse  point  of 
(3.1),  so  that  (5.13)  and  (5.14)  hold,  and  G(x,)  ^  rge  A(x,)|j.^  by  (5.6)  and  the 
nonsingularity  of  (3.1).  Then  it  follows  from  (5.16)  that  Gi((f,)  ^  rge  Ai(^,),  and  hence 
that  1^,  =  v^“^(x,)  is  a  standard  singular  point  of  the  reduction  (3.18)  of  (3.x;  near  r, 
since  Aj  in  (5.12)  and  Gi  in  (5.15)  are  exactly  as  in  that  reduction.  Conversely,  if  if,  is  a 
standard  singular  point  of  the  reduction  of  (3.1)  near  x,,  then  (5.13)  and  (5.14)  hold  (and 
also  G(x»)  ^  rge  A(x,)|^^  by  (5.16)),  and  hence  (5.6)  and  (5.7)  also  hold.  This  means 
that  X*  is  an  impasse  point  of  (3.1). 

If  X*  is  an  impasse  point  of  (3.1);  that  is,  if  =  <p~^(x»)  is  a  standard  singular  point 
of  the  reduction  (3.18)  of  (3.1)  near  x*,  it  is  straightforward  to  check  that  the  criteria  for 
accessibility  and  inaccessibility  in  Definitions  5.1  and  2.4  are  equivalent.  □ 

Remark  5.2:  From  the  above  proof  we  see  that,  by  nonsingularity  of  the  DAE  (3.1), 
the  condition  Gi(^, )  ^  rge  Ai(4,)  is  always  satisfied  when  Ai(^, )  is  singular  and  hence 
need  not  be  checked  again.  In  other  words,  to  prove  that  is  a  standard  singular  point 
of  the  reduction  of  the  nonsingular  DAE  (3.1)  (with  index  1)  near  x,,  it  suffices  to  show 
that  conditions  (2.3)  and  (2.5)  of  Definition  2.1  hold  with  ^  □ 

As  a  by-product,  Lemma  5.1  yields  an  equivalent  and  useful  definition  for  impasse  points 
of  nonsingular  DAE’s.  This  definition  is  useful,  for  it  rzirely  happens  that  intrinsic  deriva¬ 
tives  can  l)e  calculated  without  making  explicit  use  of  trivializations  and  local  coordinates. 
Definition  5.1  is  mostly  theoretical,  but  is  has  the  advantage  of  being  intrinsic:  that  is. 
independent  of  any  reduction  procedure. 

Motivated  by  Lemma  5.1  and  the  definition  of  solutions  of  singular  ODE’s  in  the  vicinity 
of  stjuidard  singular  points  (see  Definition  2.3),  we  introduce  the  concept  of  a  solution  of 
a  nonsingular  D.^E  with  index  1  near  an  impasse  point  as  follows: 
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Definition  5.2.  Let  the  DAE  (3.1)  be  nonsingular  with  index  1  and  x*  an  impasse  point 
of  (3.1 ).  For  given  6  R,  a  solution  of  the  initial  value  problem 


(5.17) 


A{x)x  =  G{x), 

x(<,)  =  X,, 


is  a  continuous  function  x  :  J  R"  on  an  interval  J  =  +  T)  or  J  =  (t*  —  T,  t*]  for 

some  T  >  0  which  is  of  class  on  =  J  \  {t,}  and  satisfies  .4(x(t))x(t)  =  G{x{t))  for 


t  G  J°. 


It  is  straightforward  to  check  that  the  reduction  procedure  of  Section  3  transforms  the 
solutions  of  (5.17)  in  the  above  sense  into  solutions  of 


\ 

in  the  sense  of  Definition  2.2.  where  Ai(^X  —  reduction  of  (3.1)  near  x*.  and 

vice-versa.  Thus,  combining  Theorem  2.1  and  Lemma  5.1,  we  find  at  once: 

Theorem  5.1.  Let  the  DAE  (3.1)  be  nonsingular  with  index  1  and  let  x,  be  an  accessible 
(resp.  inaccessible)  impasse  point  of  (3.1).  Then,  given  6  R,  the  initial  value  problem 
(5. IT)  has  exactly  two  solutions  in  the  sense  of  Definition  5.2,  both  defined  either  on 
.7  =  4-  T)  (resp.  J  =  {t,  —  T,t»])  for  some  T  >  0.  Moreover,  lim<_t,  |lx(t)||  =  oc. 

Remark  5.3:  Naturally,  the  solutions  found  in  Theorem  5.1  are  solutions  of  (5.17) 
on  J”  =  (t*,t*  -|-  T)  or  =  (t,  —  T,t^)  and  hence  are  automatically  on  this  interval 
(sec  Remark  4.1).  Moreover,  we  have  x(J°)  C  R’l  =  tt{TW  Pi  M)  by  Theorem  4.1. 
This  shows  that  impasse  points  of  (3.1)  lie  in  the  closure  of  ITi  in  IF.  In  fact,  they  lie 
on  the  boundary  of  IF]  in  IF  since  x,  €  IFj  is  impossible  for  an  impasse  point  by  the 
nonsingularity  of  (3.1).  □ 

The  behavior  of  the  solutions  of  (5.17)  in  the  vicinity  of  an  impasse  point  x,.  described 
in  Theorem  5.1,  justifies  the  terminology  “impasse  point",  at  least  in  the  acce.ssible  case, 
since  such  points  are  reached  in  finite  time  by  trajectories  emanating  elsewhere  in  IF 
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which  cannot  be  continuously  extended  beyond  that  impasse  point.  A  justification  of  the 
terminology  in  the  inaccessible  case  can  be  seen  in  the  fact  that  inaccessible  impasse  points 
become  accessible  by  simply  reversing  the  evolution  in  time. 

Because  solutions  of  a  (not  necessarily  nonsingular)  DAE  (3.1)  lie  in  iri  =  7r{W  H 
TM),  their  closure  relative  to  the  open  set  L’”  must  lie  in  W  when  W  is  closed  in  t  ",  as,  for 
example  when  the  set  (x  G  t/”  :  rank  A(x)  =  r}  is  closed  in  17"  (Proposition  3.1).  This  is 
the  case  in  many  practical  applications  but  not  the  only  mathematically  sound  possibility. 
When  W  is  not  closed  in  it  becomes  possible,  a  priori  at  least,  for  points  x,  ^  IT  to 
be  reached  in  finite  time  by  C'  trajectories  that  cannot  be  continuously  extended  beyond 
X,.  Such  points  could  also  be  called  impasse  points  but  they  are  not  covered  by  our 
theory.  Likewise,  there  may  be  points  x,  6  IT  corresponding  to  higher  singularities  of 
the  reduction;  that  is,  with  dimker  Ai(^, )  >  1,  at  which  trajectories  stop.  Examples 
of  this  kind  are  easily  obtained  by  considering  systems  of  DAE’s  in  independent  variables 
all  having  at  (say)  0  an  accessible  impasse  point  in  the  sense  of  Definition  5.1.  This  is  to 
say  that  our  treatment  of  impasse  points  of  DAE’s  is  not  exhaustive  and  that  the  impasse 
points  of  this  paper,  definitely  the  most  frequently  encountered  ones,  should  perhaps  be 
called  “standard  impasse  points”  since  other  types  may  exist  in  physical  problems. 

6.  Impasse  Points  and  Foldpoints. 

A  simple  class  of  DAE's  (3.1)  has  the  form 

f  9(3-1, X2)  =  0  G 

(6.1)  { 

[  X2  =  /(X1.X2)  G  R”. 

where  x  =  (xi,X2)  G  R"“”  x  R"  and  /  :  U”  —>■  R”,<7  :  T"  — *■  R""”  are  mappings. 
(The  differentiability  assumption  is  retained  here  for  consistency  with  the  previous  setting 
although  it  would  suffice  to  assume  only  that  /  G  C’ ). 

As  we  shall  see,  for  (6.1)  the  impa.sse  points  of  the  previous  section  coincide  essentially 
with  the  simple  foldpoints  of  ry~'(0)  relative  to  the  splitting  R"  =  R"~^  G-  R'".  This  result 
relates  our  work  to  various  discussions  of  impasse  points  found  in  the  literature  where  the 
form  (6.1)  is  assumed,  notably  [CDS9]  (see  also  [M91]). 
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Theorem  6.1.  ('j  The  DAE  (6.1 )  satisfies  the  general  hypothesis  of  Section  3  if  and  only 
if  the  mapping  g  is  a  submersion  in  ?7". 

(ii)  If  (i)  holds  then  (6.1)  is  nonsingular  with  index  1  if  and  only  if 

(6.2)  {x  e  g~\0),  D2g(x)f{x)  e  rge  Dig(x)}  Dig{x)  e  GLiW~^), 

where  Da  denotes  partial  differentiation  with  respect  io  Xq,  a  =  1,2. 

(Hi)  If  both  (i)  and  (ii)  hold  then  x*  G  U"  is  an  impasse  point  of  (6.1)  if  and  only  if 
X,  G  and 

(6.3)  dimkerZ)i5{x,)  =  1, 

(6.4)  (D?g(x.)(u)^«)  ^0, 

for  some,  or  equivalently  all,  pairs  of  nonzero  vectors  (u,u)  in  ker  Digf(x, )  x  ker  D\g{x^)^ , 
where  {■,  •)  denotes  the  natural  inner  product  of  R”"'".  Moreover,  the  impasse  point  x,  is 
accessible  or  inaccessible  if  and  only  if 

(6.5)  a(x,)(u,{j)  =  {D2g(x„)f(x»),u){D'^g{x»){u,u),u) 

is  <  0  or  >  0,  respectively,  for  some,  or  equivalently  every,  pair  (u,u)  of  nonzero  vectors 
in  ker  Dig(xt)  x  ker  Hi  ^(x.  )^. 

Proof:  The  DAE  (6.1)  is  the  special  case  of  (3.1)  in  which 

/O  0\ 

(6.6)  .4(x)  =  I  I  ,  H  =  identity  of  R*", 

VO  IrJ 


is  independent  of  x  and 


(6.7) 


G  R"-’'  X  R^ 
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For  the  proof  of  (i)  we  write 


P  = 


F{x,p)  =  A{x)p  -  G{x)  = 


Then  for  (/i,  5)  €  R"  x  K."  we  have 


DF{x,p){h,q) 


/  Dg{x)  \ 
\q2  -  Df{x)h  ) 


which  makes  it  obvious  that  F  is  a  submersion  on  x  R"  if  and  only  if  5  is  a  submersion 
on  L’".  In  this  case,  M  =  F~^{0)  =  {(ar,p)  6  x  R";  x  £  g~^i0),p2  =  fi^)}  has 

the  projection  W  =  ■n[M)  =  Since  A{x)  is  independent  of  x  and  has  rank  r  the 

condition  rank  A{x)  =  r,  Va*  G  T'F,  holds  trivially. 

For  (ii)  assume  that  5  is  a  submersion  on  f/”.  Because  of  W  =  ^“^(0)  we  have 


(G.S) 


a:  G  IF  =>  T^W{h  = 


G  R”  X  R'' :  Dig{x)hi  4-  D2g{x)h2  =  0} 


and  from  (6.6)  and  (6.8)  we  find  that 


(6.9) 


ker>l(a:)|j.^j^  =  kerZ?i(/(x)  x  {0} 


and 


(6.10)  rge  ^(x)|^^^,.  =  { 


^  G  R"""  X  R"  :  3hi  G  R'^-\Dig{x)hi  +  D2g(x)h2  =  0}. 


Thus,  it  follows  from  (6.7)  and  (6.10)  that  x  G  IF  =  fir~^(0)  and  G{x)  G  rge  4(x)|^^^^.  if 
and  only  if  x  G  g~^(0}  and  D2g(x)f{x)  G  rge  D\f{x).  Also,  since  T^W  is  r-dimensional, 
rank  -4(x)|.^^^^,  =  r  if  and  only  if  A{x)\^^^^.  is  one-to-one;  that  is,  by  (6.9)  if  ker  Dig{x)  = 
{0}.  In  turn,  this  amounts  to  D\g{x)  G  GL(R"~’'),  and  hence  the  condition  (6.2)  is  indeed 
equivalent  to  the  DAE  (6.1)  being  nonsingular  with  index  1  (see  Definition  4.1). 
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Finally,  for  the  proof  of  (iii)  let  x*  6  U"  be  an  impasse  point  of  (6.1).  Then,  by 
definition,  x,  G  IF  =  gr“^(0).  To  prove  that  conditions  (6.3)  and  (6.4)  characterize  x,  as 
an  impasse  point,  we  use  Lemma  5.1  to  show  that  (6.3)  and  (6.4)  are  equivalent  to  the 
condition  that  if*  =  (^~^(x*)  is  a  standard  singular  point  of  the  reduction  of  (6.1)  near  x, 
(under  the  standing  assumption  that  (6.1)  is  nonsingular  with  index  1). 

With  if  I  :  V  W  =  5f“^(0)  being  a  diffeomorphism  of  an  open  subset  U’^  of  R'' 
onto  an  open  neighborhood  O  of  x,  in  IF,  the  method  of  Section  3  with  A  and  G  given 
by  (6.6)  and  (6.7),  respectively,  shows  that  the  reduction  of  (6.1)  near  x,  has  the  form 

(6.11)  = 

with 

(6.12)  ,4i(f)  =  P.T>9(f) 

where  P,  is  the  projection  onto  R'"  corresponding  to  the  splitting  R"  =  R"~’'  x  R*”  and 

(0.13)  GAO  =  p.G{v(m- 

With  f,  =  (^"^(x,)  we  have  €  Gi(R^, Ti. PF).  Thus,  b}'  (6.8)  and  (6.12)  it 

folio  .vs  that 

(6.14)  ker-4i(f,)  =  Pv^(f,)“’ (ker  Pi5^(x.)) 

and  therefore  that  dimker,4i(f, )  =  1  if  and  only  if  dim ker  Z)i<7(x* )  =  1;  that  is,  if  (6.3) 
holds.  Now,  we  show  that  (6.4)  is  equivalent  to  the  condition 

(6.15)  (P,4i(f*)7/)7?  rge.4i(f.), 

where  j]  G  ker.4i(f,)  \  {0}  and  thus,  by  (6.14),  r]  =  £)(p(f,)“^u,M  G  kev  D\g(Xt,)  \  {0}. 
Bccau.se  of  dim  ker  ,4 1  ( f » )  =  1,  (6.15)  amounts  to 

(6.16)  ((P>.4i(f.)7/)77  ,77)  ^0, 
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where  ^  is  any  nonzero  vector  in  (rge  Ai(^* ))■*■.  We  claim  that  any  such  vector  has 
the  form  fj  —  D2g{x*)^u  with  u  G  kerDi$f(x*)  \  {0}.  To  see  this,  note  first  that 
Dg{x^)^  G  is  one-to-one  since  Dg{x»)  G  £(R",R"“'')  is  onto.  Thus,  for 

u  G  kerDip'(x*)^  \  {0},  we  have  D2g{x^.)^u  ^  0.  Next,  from  (6.8)  and  )  G 

GZ(R’',  Tj.^  T'T)  we  deduce  that 

(6.17)  rge  Ai(^.)  =  {/i2  G  R*"  :  3hi  G  R"'"', T>iff(a-, )/ii  +  £>2fif(x,)/i2  =  0}. 

It  follows  at  once  that  D2g{xt)^u  G  (rge  A]((f,))‘‘"  for  all  u  G  ker  Dig{x^)^  \  {0}.  This 
proves  the  claim  since,  evidently,  D2g(x^)^u  ^  0  and  (rge  Ai(if, ))-*-  is  one-dimensional. 
From  (6.12)  we  deduce  that  (6.16)  may  be  written  as 

(6.18)  {P,D'^(p{^^){r],T]),(D2gix»))'^u)  7^  0, 

for  T]  =  and  (u,u)  an  arbitrary  pair  of  nonzero  vectors  in  kerDi^(x,)  x 

ker  Di<7(x,  )^.  To  calculate  the  left-hand  side  of  (6.18)  we  use  the  relation  goip  =  0  follow¬ 
ing  from  the  definition  of  (,?.  Differentiating  twice  at  ^  we  obtain  D^g{x»){D^{^»)ri)'^  + 
Dg{x»)D'^ip(x»)(7])^  =  0.  Taking  now  the  inner  product  with  u  and  using  g  = 
we  get  (I?'(,?(^. )(//,  7?),  (Z)(/(x, ))^i/)  =  (D^g(x*)(u,u),u).  But  the  left-hand  side  is  just 
(P,Dy((f,  1(77,  r?)  {D2g{x^))^u)  because  {D\g(xt,))^u  —  0  and  the  right-hand  side  is 
—  {D\g{x^)(u.v),  u)  because  cf  u  G  kerDig(x,)  C  R""^.  Thus, 

(6.19)  {P,D-^p{^„)(g.g).D2g{x^)'^u)  =  -{D‘\g{x»){u,u),u) , 

which  shows  that  (6.4)  is  equivalent  to  (6.18)  and  hence  to  (6.15). 

By  Remark  5.2  the  equivalence  of  (6.3)  with  dimker  A](^,)  =  1  and  of  (6.4)  with  (6.15) 
suffices  to  prove  that,  together,  (6.3)  and  (6.4)  are  equivalent  to  x,  being  an  impasse  point 
of  (6.1). 

For  th('  accessibility  or  inaccessibility  criterion  note  that  by  the  previous  calculation 
{[DA\(S,*)g)g.ri)  =  —{D^g{x»)(u,u),u)  and  that  Iw  (6.18)  and  -fj  =  D2g{x,)^u  we  have 
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(Gi(i^*),77)  =  {fixt)-,D2gixt)^u).  From  these  relations  and  the  criterion  (2.6)  for  accessi¬ 
bility  or  inaccessibility  of  we  obtain  the  criterion  (6.5)  for  the  accessibility  or  inaccessi¬ 
bility  criterion  of  x*.  □ 

Under  the  sole  assumption  that  </  is  a  submersion  in  U”  the  conditions  (6.3)  and  (6.4) 
characterize  x*  G  ^~^(0)  aa  a  simple  foldpoint  of  ^“^(0)  relative  to  the  splitting  R”  — 
R""""  X  R'’;  that  is,  as  a  turning  point  in  the  case  when  r  =  n  —  1  and  hence  g  is  scalar.  The 
accessibility  or  inaccessibility  criterion  (6.5)  has  no  geometric  interpretation  and  translates 
a  purely  dynamic  property. 

In  [CDS9],  Chua  and  Deng  characterize  an  impasse  point  of  (6.1)  as  a  point  x,  = 
{x»i,Xt2)  such  that  the  mapping  (x^.A)  i — >  <7(a:i,x,2  +  has  a  turning  point  at 

('^i*,0).  Now,  a  simple  (or  generic)  turning  point  (xi,,0)  of  this  mapping  is  charac¬ 
terized  by  the  three  conditions:  dimker Dig(x,)  =  1,  D2g(x^)f{x,)  ^  rge  Dig{x^)  and 
{D\g{x»  ){u)'^ ,u)  ^  0  for  every  pair  of  nonzero  vectors  (u,u)  6  kerZ)ig(x,  )  x  ker  Dig{x»)^. 
In  our  exposition,  the  first  and  third  of  these  conditions  are  exactly  (6.3)  and  (6.4),  re¬ 
spectively,  and  the  second  one  is  contained  in  (6.2)  and  (6.3).  Thus,  for  the  special  case 
of  (6.1)  our  definition  of  impasse  points  is  essentially  the  same  as  that  of  [CD89]. 

7.  Examples. 

7.1  A  nonlinear  LPC  network.  Consider  a  nonlinear  LRC  network,  consisting  of  a 
resistor  branch  (1),  a  capacitor  branch  (2)  and  an  inductor  branch  (3)  in  parallel.  Denote 
by  I,  the  currents  in  the  U*-branch,  and  by  V  the  (common)  voltage  drop  in  each  branch. 
Let  C(U)  and  Lils)  represent  the  capacity  mid  the  inductivity,  respectively,  where  C  : 
R  — >■  (0,  oc)  and  L  :  R  — >  (O.  oo)  are  smooth  enough  functions  (see  e.g.  Smale  [S72]  or 
Takeiis  [T7C]).  From  Kirchhoff’s  law,  we  have  I\  +  I2  +  I3  —  0.  while  the  time-evolution 
of  the  voltage  drop  and  currents  in  the  capacitor  and  inductor  branches  is  given  by  I2  = 
C{V)V .  V'  =  T(/3)/3.  Finally,  in  the  resistor  branch  (1),  /]  and  V  are  related  through  an 
equation  of  the  form  g(I\ .  U)  =  0,  where  (/  :  R  x  R  — +  R  is  a  smooth  enough  function. 

For  convenience,  we  shall  use  the  variables  .x,  ~  Ii,  1  <  ?  <  3,  x.i  =  U.  In  this  notation, 
the  (xjuations  governing  the  evolution  of  the  voltage  drop  and  currents  of  the  above  network 
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have  the  form 


(7.1) 


-4(3’)ir  =  G{x) 


where  x  —  (a:i ,  a’2,  X3,  X4  and 


i .. 


-4(x) 


/O  0  0  0  \ 

0  0  0  0 

0  0  L(x3)  0 

Vo  0  0  C(X4)/ 


,  G(x) 


/  X\  4-  X2  -f  X3  \ 

?(Xl,X4) 

X4 

V  X2  / 


Clearly,  rank  A{x)  =  2,  Vx  €  since  L  and  C  assume  strictly  positive  values  by 
hypothesis.  Under  the  mild  condition  that  (dg/dx4){xi,  X4)  ^  0  for  x  G  R"*,  the  reduction 
procedure  of  Section  3  can  be  applied.  In  fact,  here  the  manifold  TU  is  given  by  W  =  {x  G 
R'*  :  X]  +  X2  +  X3  =  0,  5'(xi,X4)  =  0},  and  hence  can  be  parametrized  by  if  =  (X],X2), 
i.e.  ir  =  {(xi,X2,— xi  —  X2,0(xi))  :  (xi,X2)  €  R^},  where  ^  :  R  R  is  characterized  by 
_(7(xi .  ^(xi ))  =  0.  We  leave  it  to  the  reader  to  check  that  the  DAE  (7.1)  is  nonsingular 
(Definition  4.1)  in  some  open  subset  U*  C  R^  if  and  only  if 


(7.3) 


{xG  iunu\— (xi) -0}  =>X2  ^0. 
dx  1 


and  tliat  the  reduced  system  reads 


(7.4) 


Ai(Xi,X2)  j  )  =  Gi(xi.X2), 
X2 


vlierc' 


.4|(.r].X2)  = 


L(-xi  -  X2)  L{-xx  -  X2) 


0 


C(^(x,))i^(.r,) 


,  G(x,.X2) 


-(9(xi ) 
.72 


aiK 


From  Lemma  5.1.  x  =  {x\,  X2.  —x\  —  X2,d(x{))  G  IF  Cl  I'  is  an  impasse  jjoint  of  ( 7.1 )  if 
1  only  if  (x  1 ,  X2  )  is  a  standard  singular  point  of  (7.4):  that  is.  -^{x\  )  =  0  .  ■f77(.?'i )  ^  0. 
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Under  the  assumption  (7.3)  this  is  straightforward  to  verify;  see  also  Remark  5.2.  Moreover, 
accessibility  or  inaccessibility  of  the  impasse  point  depends  upon  <  0  or  >  0, 

respectively.  For  instance  if  g(xi,X4)  =  X4  —  Xj  —  y,  y  ^  R  (an  example  chosen  in 
Takens  [T76]),  then  ^(xi)  =  7  +  Xj  and  the  impasse  points  of  (7.1)  axe  points  of  the  form 
(0, X2,  — X2, 7),  X2  0  (accessible  if  X2  <  0,  inaccessible  if  X2  >  0).  The  DAE  (7.1)  is 
singular  in  any  open  subset  U‘^  containing  the  point  (0,0, 0,7),  for  (7.2)  fails  to  hold  and, 
in  fact,  this  point  is  not  an  impasse  point  but  a  “fimnel”  if  7  ^  0  (see  [T76];  see  also 
[RR92]  for  further  comments  and  numerical  results). 

Remark  7.1:  The  choice  of  A{x)  and  G(x)  in  (7.2)  is  not  canonical.  For  instance,  (7.1) 
is  unchanged  if,  instead  of  (7.2),  we  choose 


A(x) 


/O  0  0  0\ 
0  0  0  0 
0  0  10 
Vo  0  0  1/ 


,  G(x) 


/Xi  +  X2  +  X3\ 

g(xi,X4) 

X4/L{Xi) 

\  X2/C(X4)  / 


With  the  above  choice,  the  problem  fits  into  the  framework  of  Section  6  and  hence  the 
impasse  points  of  (7.1)  appear  as  foldpoints  of  the  mapping 


(Xj,X2,X3,X4) 


(a-1  +  X2  +  X3 
g{xi , X4 ) 


relative  to  the  splitting  =  R?  ,  x  R^  ,.  □ 

7.2  One-parameter  stationary  problems.  We  shall  show  that  one-parameter  sta¬ 
tionary  problems,  in  which  no  evolution  of  any  kind  is  involved,  can  be  reformulated  as 
a  quasilinear  DAE,  usually  of  index  <  1,  and  that  their  turning  points  coincide  with  the 
impasse  points  of  this  paper.  In  fact,  all  the  singularities,  bifurcations  and  others,  that 
may  appear  in  a  stationary  one-parameter  problem  may  also  be  characterized  as  singular¬ 
ities  of  the  DAE  formulation.  Thus,  the  analysis  of  DAE's,  their  impasse  points  and  their 
singularities  encompasses  the  entire  study  of  on^parameter  stationary  problems,  their 
turning  points  and  singularities. 
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Let  us  begin  with  the  equation 


(7.5)  /(A,x)  =  0, 

where  /  :  R  X  R"  — +  R"  is  a  smooth  enough  function.  If  solutions  are  sought  in  the  form 
X  =  x(A),  a  function  of  the  parameter  A  (which  is  legitimate  as  long  as  the  derivative 
Dif{X,x)  is  invertible)  then  x(A)  is  a  solution  of  Dxf{X^x)  +  Dzf(X,x)x  =  0,  where 
the  dot  represents  d/dX.  Introducing  the  auxiliarj'  equation  A  =  1,  this  equation  can  be 
rewritten  as 

where 

A(A,x)=(^  ^  ° 

\Dxf(X,x)  Dxf(X,x)J 

Evidently,  .4(A,x)  is  singular  if  and  only  if  Dif{X,x)  is  singular,  and  a  straightforward 
verification  shows  that  {X,x)  is  a  standard  singular  point  of  (7.6)  if  and  only  if 

(i)  dimker  Di/(A,  r)  =  1, 

(ii)  Dxf{X,x)^  rgeDi/(A,r), 

(iii)  Dlf{X,x)h'^  ^  rge  Dxf{X,x),  'ih  eker  Dxf{X,x)  -  {0}. 

As  is  well-known,  conditions  (i)  -  (iii)  above  characterize  the  solution  (r,A)  of  (7.5)  as  a 
simple  turning  point  of  (7.5).  Since  if  /  is  constant  along  the  trajectories  of  (7.5),  standard 
singular  points  of  (7.6)  along  trajectories  emanating  at  (Ao,xo)  with  /(Ao,xo)  =  0  coincide 
with  the  simple  turning  points  of  (7.5).  More  generally,  equivalence  between  singularity 
of  A{X,x)  and  singularity  of  Dj/(A,r)  shows  that  all  the  singularities  of  (7.5)  appear  as 
singularities  of  the  quasilinear  singular  ODE  (7.6). 

Quasilinear  DAE’s  rather  than  ODE’s  are  needed  to  reformulate  constrained  problems, 
i.e.  problems  of  the  form 

[  (j{x)  =  0  G  R"“'', 

(7.7)  { 

f(X,x)  =  0  €  RL 
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In  (7.7),  the  equation  g{x)  —  0  may  be  viewed  as  a  condition  characterizing  the  admissible 
“state”  variables  x.  Using  the  same  procedure  as  before,  i.e.  differentiating  the  identity 
/(A,  2:(A))  =  0,  we  see  that  (7.7)  can  be  reformulated  as  the  problem  of  solving  the  DAE 

(7.8)  A(A,x)  ^  ^  =  G(A,x), 

where  .4(A,  x)  is  the  (n  +  1)  x  (n  +  1)  matrix 

.1  0  . 

A{X,x)=  DxfiKx)  D,f{Kx) 

\  0  0  / 

and  G(A,x)  €  R  x  R”  x  R”"’'  =  is  given  by  G(A,x)  =  (1,0, ^(x))^  (independent 

of  A). 

It  is  easily  seen  that  under  the  general  conditions  ensuring  that  (7.8)  is  a  nonsingular 
DAE,  the  impasse  points  of  (7.8)  along  trajectories  emanating  at  a  point  (Ao,xo)  where 
g(xo)  =  0,  and  /(Ao,xo)  =  0,  coincide  with  the  simple  turning  points  of  the  equation 
—  0,  where  ;  R”  — >  R”  is  a  local  parametrization  of  the  manifold  $'“^0)  (Note 
that  5'~^(0)  is  a  manifold  because  (7.8)  is  a  nonsingular  DAE). 

8.  Higher  Index  Problems. 

It  was  shown  in  Section  3  that  under  Assumption  3.1,  the  n-dimensional  quasilinear 
DAE  in  7/"  C  R" 

(8.1 )  A(x)x  =  G(x) 
can  be  reduced  locally  to  the  r-dimensional  form 

(5.2)  Aj(a^  =  Gi(a 

where  r  <  n  is  the  (constant)  rank  of  A(x).  Recall  that  (8.1)  was  called  a  nonsingular 
DAE  with  index  1  if  the  conditions  of  Definition  4.1  are  fulfilled.  If  so,  (8.2)  reduces  to 
the  ODE  =  .4i((^)“’ Gi(i^)  (see  the  proof  of  Theorem  4.1). 
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Now,  it  may  also  happen  that  (8.2)  is  a  nonsingular  DAE  with  index  1.  In  this  case. 
(8.1)  will  be  called  a  nonsingular  DAE  with  index  2.  From  the  reduction  procedure  of 
Section  3,  it  appears  that  if  (f,  is  an  impasse  point  of  (8.2),  the  images  x(t)  =  of 

the  two  solutions  of  (8.2)  emanating  or  terminating  at  (f,  (in  the  notation  of  Section  3) 
solve  (8.1)  and  both  emanate  or  terminate  at  ar,  =  Also,  no  other  solution  of  (8.1) 

passing  through  a:,  exists.  Thus,  x,  possesses  all  the  properties  characterizing  impasse 
points  of  quasilinear  DAE’s  of  index  1  and  hence  should  be  referred  to  as  an  impasse  point 
of  the  quasilinear  DAE  (8.1)  of  index  2.  It  is  straightforward  (but  somewhat  cumbersome) 
to  check  that,  as  expected,  this  definition  is  independent  of  the  specific  parametrization  y; 
used  to  reduce  (8.1)  to  the  form  (8.2).  It  is  equally  clear  how  higher  index  (quasilinear) 
DAE’s  and  their  impasse  points  can  be  defined.  Naturally,  these  definitions  implicitly 
require  sufficient  smoothness  of  A  and  G  in  (8.1).  Note  that  even  if  -4(x)  in  (8.1)  is 
independent  of  x,  Ai(i^)  in  (8.2)  will  not  be  constant  in  general.  Thus,  for  higher  index 
problems,  the  simplified  framework  of  Section  6  cannot  be  used. 
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